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Abstract 

We evaluate differential distributions for exclusive scalar /o(1500) meson (glueball candidate) 
production for pp iViiV2/o (FAIR@GSI) and pp ppfo (J-PARC@Tokai). Both QCD diffrac- 
tive, pion-pion meson exchange current (MEC) components as well as double-diffractive mechanism 
with intermediate pionic loop are calculated for the first time in the literature. The pion-pion com- 
ponent, which can be reliably calculated, dominates close to the threshold while the diffractive 
component may take over only for larger energies. At the moment only upper limit for the QCD- 
diffractive component can be obtained. The diffractive component is calculated based on two-gluon 
impact factors as well as in the framework of Khoze-Martin-Ryskin approach proposed for diffrac- 
tive Higgs boson production. Different unintegrated gluon distribution functions (UGDFs) from 
the literature are used. Rather large cross sections due to pion-pion fusion are predicted for PANDA 
energies, where the gluonic mechanism is shown to be negligible. The production of /o(1500) close 
to threshold could limit the so-called ttNN form factor in the region of larger pion virtualities. We 
discuss in detail the two-pion background to the production of the /o(1500) meson. 
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I. INTRODUCTION 



Many theoretical calculations, including lattice QCD, predicted existence of glueballs 
(particles dominantly made of gluons) with masses M > 1.5 GeV. No one of them was 
up to now unambiguously identified. The nature of scalar mesons below 2 GeV is also 
not well understood. Lattice QCD approach with quenched quarks find a scalar gluonium 
(glueball) at approximately 1.6 GeV [l| ^. Also the analyses in the framework of chiral 
LagrangianslsL 4 1 indicate that /o(1500) is dominantly gluonium state. The QCD Sum 
Rules p, 0, oris, a, l3| suggest that the states at approximately 1 GeV and 1.5-1.6 GeV 



are admixtures of gluonium and qq states. A recent analysis in the framework of Gaussian 
QCD Sum Rules [llj, which is well suited to qq - gluonium mixing, find that the states at 
about 1 GeV (/o(980)) and at about 1.4 GeV are strongly mixed with the preference of the 
higher-mass state to have slightly larger gluonium admixture. Summarizing this discussion, 
it may be very difficult to find a clear signal of gluonium. Further studies of the scalar meson 
production in several processes may shed more light on the quite complicated problem. 

The lowest mass meson considered as a glueball candidate is a scalar /n(1500 ) |12|] dis- 
covered by the Crystall Barrel Collaboration in proton-antiproton annihilation [l3|. The 
branching fractions are consistent with the dominant glueball component [13] • It was next 
observed by the WA102 Collaboration in central production in proton-proton collisions in 
two-pion [15I] and four-pion [l6| decay channels at \fs ~ 30 GeV ^. Close and Kirk [iT^ pro- 
posed a phenomenological model of central exclusive /o(1500) production. In their language 
the pomerons (transverse and longitudinal) are the effective (phenomenological) degrees of 
freedom [l8|. The Close-Kirk amplitude was parameterized as ^ 



A^(ti, t2, </>')= AT exp (^y (ti + ts)^ +aL^^exp(^y(ti+t2)j cos((/.'). (1.1) 

In their approach there is no explicit /o(1500)-rapidity dependence of the corresponding 
amplitude. Since the parameters were rather fitted to the not-normalized WA102 data 
[l5| no absolute normalization can be obtained within this approach. Furthermore the 
parameterization is not giving energy dependence of the cross section, so predictions for 
other (not-measured) energies are not possible. In the present paper we will investigate 
rather a QCD-inspired approach. It provides absolute normalization ^, energy dependence 
and dependence on meson rapidity (or equivalently on of the meson). 

The nature of the /o(1500) meson still remains rather unclear. New large-scale devices 
being completed (J-PARC at Tokai) or planned in the future (FAIR at GSI) may open a 
new possibility to study the production of /o(1500) in more details. 

In the present analysis we shall concentrate on exclusive production of scalar /o(1500) in 



^ The approaches with dynamical quarks find relatively large mixing with qq states 

^ No absolute normalization of the corresponding experimental cross section was available. Only two-pion 

or four-pion invariant mass spectra were discussed. 
^ The 4>' dependence applies in the meson rest frame (current-current cm.) 
^ As will be discussed later it is rather upper limit which can be easily obtained. 
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the following reactions: 



p + p^p + foil500)+p , 

P + P^P + fo{1500)+p, (1.2) 
p + p n + fo{1500) + n . 

While the first process can be measured at J-PARC, the latter two reactions could be mea- 
sured by the PANDA Collaboration at the new complex FAIR planned in GSI Darmstadt. 
The combination of these processes could shed more light on the mechanism of /o(1500) 
production as well as on its nature. 

If /o(1500) is a glueball (or has a strong glueball component p^]) then the mechanism 
shown in Fig. [1] may be important, at least in the high-energy regime. This mechanism 
is often considered as the dominant mechanism of exclusive Higgs boson 2^ and Xc{0~^) 



meson 



21| production at high energies. There is a hope to measure these processes at LHC 
in some future when forward detectors will be completed. At intermediate energies the same 
mechanism is, however, not able to explain large cross section for exclusive rj' production 
2^ as measured by the WA102 Collaboration. Explanation of this fact is not clear to us in 
the moment. 




FIG. 1: The sketch of the bare QCD mechanism for diffractive production of the glueball. The 
kinematical variables are shown in addition. 



At lower energies {^/s < 20 GeV) other processes may become important as well. Since 



the two-pion channel is one of the dominant decay channels of /o(1500) (34.9 ± 2.3 %) ^ 
one may expect the two-pion fusion (see Figj2]) to be one of the dominant mechanisms of 
exclusive /o(1500) production at the FAIR energies. The two-pion fusion can be also relative 
reliably calculated in the framework of meson exchange theory. The pion coupling to the 
nucleon is well known The ttNN form factor for larger pion virtualities is somewhat less 
known. This may limit our predictions close to the threshold, where rather large virtualities 
are involved due to specific kinematics. At largest HESR (antiproton ring) energy, as will be 
discussed in the present paper, this is no longer a limiting factor as average pion virtualities 
are rather small. 

In this paper we concentrate on the mechanism of the reaction. Our aim here is to 
explore a possibihty of studying exclusive /o(1500) meson production in the FAIR and J- 
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FIG. 2: The sketch of the pion-pion MEG mechanism. Form factors appearing in different vertices 
and kinematical variables are shown exphcitly. 




FIG. 3: The sketch of the double-diffractive mechanism with pionic loop for exclusive production 
of the glueball candidate /o(1500). The stars attached to vr mesons denote the fact that they are 
off-mass-shell. 



PARC energy range and explore the potential of these facilities. While their are some ideas 
about the reaction mechanism at higher energies, the mechanism at lower energies was never 
studied. We shall investigate new mechanisms of pion-pion fusion shown in FigJJl the QCD 
mechanism shown in Fig{T] and a mechanism with intermediate pionic loop shown in FigJ31 
The second (QCD) mechanism is typical for high energies but here we wish to investigate its 
role at intermediate energies and in particular its vanishing at low energies and the interplay 
with the pion-pion fusion mechanism. 
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II. EXCLUSIVE PROCESSES 



A. Cross section and phase space 

The cross section for a general 3-body reaction pp pp/o(1500) can be written as 

c^^„M = W ■ d'PS . (2.1) 

Above m is the mass of the nucleon. 

The three-body phase space volume element reads 

d^PS = ■ (27r)^54(pi +P2- p[ - p', - Pm) ■ (2.2) 

2E[{2Txf2E'^{2'Kf2EM{2'KY 

At high energies and small momentum transfers the phase space volume element can be 



written as 25 



d^PS ^ ^dtidt2d^id^2d(j) S {s{l - ei)(l - 6) - M^) , (2.3) 

where ^i, ^2 are longitudinal momentum fractions carried by outgoing protons with respect to 
their parent protons and the relative angle between outgoing protons G (0, 27r). Changing 
variables (^1,^2) {xp,M'^) one gets 

d^PS ^ -^dhdt2 — '^^^ d(t) . (2.4) 

2'vr4 %^x| + 4(M2 + |Pm,P)/s ^ ' 

The high-energy formulas (12. 3p and (12.41) break close to the meson production thresh- 
old. Then exact phase space formula (12.21) must be taken and another choice of variables 
is more appropriate. We choose transverse momenta of the outgoing nucleons {Pit,P2t)j az- 
imuthal angle between outgoing nucleons (0) and rapidity of the meson (y) as independent 
kinematically complete variables. Then the cross section can be calculated as: 



d(y = 2_^J {pit,P2t,(l),y)\k ^ . ON rc,„\^nTP'nTP' nP^iP^tdpitdp2td(l)dy, (2.5) 

k 



2^s(s-4m2) (27r)5 2E[ 2E'^ 2 



where k denotes symbolically discrete solutions of the set of equations for p\^ and p'2 

Vs- Em = + p'l^ + y^^4T^ 

-Pmz = V'lz + V'2. 

where mu and m2t are transverse masses of outgoing nucleons. The solutions of 
Eg. (12.61) depend on the values of integration variables: p'^^ = Piz{PityP'2ty 'P^v) aiid = 
P2z{PuyP2ty ^^y)- The extra Jacobian reads: 



' (2.6) 

z ) 



P'lzik) p'2z{k) 



v^it+p'iziky v^it+P2ziky 



(2.7) 



In the limit of high energies and central production, i.e. p[^ ^ (very forward nucleonl), 
—p'2z ^ (very backward nucleon2) the Jacobian becomes a constant J' ^ ^. 

The matrix element depends on the process and is a function of kinematical variables. 
The mechanism of the exclusive production of /o(1500) close to the threshold is not known. 
We shall address this issue here. Therefore different mechanisms will be considered and the 
corresponding cross sections will be calculated. 
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B. DifFractive QCD amplitude 



According to Khoze-Martin-Ryskin approach (KMR) [20|, the amphtude of exclusive 
double diffractive colour singlet production pp j9p/o(1500) can be written as 

The normalization of this amplitude differs from the KMR one j2Q| by the factor s/2 and 
coincides with the normalization in our previous work on exclusive r^'-production |22|. The 
amplitude is averaged over the colour indices and over two transverse polarisations of the 
incoming gluons 20|. The bare amplitude above is subjected to absorption corrections which 
depend on collision energy (the bigger the energy, the bigger the absorption corrections). 
We shall discuss this issue shortly when presenting our results. 
The vertex factor = 

ti Plit) expression (12.81) describes the coupling of 
two virtual gluons to /o(1500) meson. Recently the vertex was obtained for off-shell values 
of qi^t and g2,t in the case of Xc(0) exclusive production \1V\. An almost alternative way 
to describe the vertex is to express it via partial decay width T{M — > gg). ^ The latter 
(approximate) method can be used also for the /o(1500) meson production. 

In the original Khoze-Martin-Ryskin (KMR) approach j20| the amplitude is written as 



M = N I ^^go,tP[/o(1500)] ^^^^^^^^ ^ t,)/f ^^^(x„ 4, Ql,^ t,) , (2.9) 



rf^g o,tP[/o(150 0)] 

where only one transverse momentum is taken into account somewhat arbitrarily as 

Q\,t = min{go,t' 9i,t} ' <52,t = min{g2^, ql ^} , (2.10) 

and the normalization factor N can be written in terms of the /o(1500) — > gg decay width 
(see below). 

In the KMR approach the large meson mass approximation M ^ |qi,t|, |q2,t| is adopted, 
so the gluon virtualities are neglected in the vertex factor 

P[/o(1500)] ~ {q,,tq2,t) = {qo,t + p[,t){Qo,t - P2,t)- (2.11) 

The KMR UGDFs are written in the factorized form: 

/™(x, x', Ql t) = f^^'^ix, x\ Ql /i^) exp(M) (2.12) 



with bo = 2 GeV~^ |20|. In our approach we use somewhat different parameterization of the 
t-dependent isoscalar form factors. 

Please note that the KMR and our (general) skewed UGDFs have different number of 
arguments. In the KMR approach there is only one effective gluon transverse momentum (see 
Eq.f l2.10l) ) compared to two independent transverse momenta in general case (see Eg. (12. 161) ). 



^ The last value is not so well known. We shall take r(M — > gg) = ^*m- This will give us an upper 
estimate. As a consequence this will allow us to show that the gluonic component is negligible for future 
experiments with the PANDA detector. 
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The KMR skewed distributions are given in terms of conventional integrated densities g 
and the so-called Sudakov form factor T as follows: 



Rr. 



d 



d In Qf 



TiQlfi^)xgix,Ql 



(2.13) 



The square root here was taken using arguments that only survival probability for hard 
gluons is relevant. It is not so-obvious if this approximation is reliable for light meson 
production. The factor Rg in the KMR approach approximately accounts for the single 
logQ^ skewed effect j20j. Please note also that in contrast to our approach the skewed 
KMR UGDF does not explicitly depend on x' (assuming x' <^ x ^ 1). Usually this 
factor is estimated to be 1.3-1.5. In our evaluations here we take it to be equal 1 to avoid 
further uncertainties. Following now the KMR notations we write the total amplitude (12. 8p 
(averaged over colour and polarisation states of incoming gluons) in the limit M ^ gi g2,i 
as 



M = Att^ 



off. 



c?'go,ti"[/o(1500)] 



•Jo,* Q'it 



(2.14) 



where the normalization constant is 



A' 



647rr(/o(1500) ^ gg) 
{Nl - \)Ml 



(2.15) 



In addition to the standard KMR approach we could use other off-diagonal distributions 
(for details and a discussion see |2l|, l22|). In the present work we shall use a few sets of 
unintegrated gluon distributions which aim at the description of phenomena where small 
gluon transverse momenta are involved. Some details concerning the distributions can be 
found in Ref. 2^. We shall follow the notation there. 

In the general case we do not know off-diagonal UGDFs very well. In 21, l^] we have 
proposed a prescription how to calculate the off-diagonal UGDFs: 



Fl{t2) 



f, 



off 
9,2 



(2.16) 



where -Fi(ti) and -^1(^2) are isoscalar nucleon form factors. They can be parameterized as 

(H) 



Fl(tl,2 



4m^-2.79ti,2 



(4m2-ti,2)(l-ti,2/0.71)2 • 



(2.17) 



In the following for brevity we shall use notation ti^2 which means ti or t2- Above ti and t2 
are total four-momentum transfers in the first and second proton line, respectively. While 
in the emission line the choice of the scale is rather natural, there is no so-clear situation for 
the second screening-gluon exchange (22| . 

Even at intermediate energies {W = 10-50 GeV) typical relatively small 

(~ 0.01). However, characteristic xi,X2 ~ Mf^^/y/s are not too small (typically > 10~^). 
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Therefore here we cannot use the small-x models of UGDFs. In the latter case a Gaussian 
smearing of the collinear distribution seems a reasonable solution: 

kl = x/°"(a;, fil) ■ Fcaussikl, ao) , (2.18) 

where (y'™"(x, /i|,) are standard collinear (integrated) gluon distribution and fcaussikf] ctq) is 
a Gaussian two-dimensional function 

Fcaussikl do) = ^ exp {-k^/2a^o) /vr . (2.19) 

Above (To is a free parameter which one can expect to be of the order of 1 GeV. Based on 
our experience in [22| we expect strong sensitivity to the actual value of the parameter cto. 
Summarizing, a following prescription for the off-diagonal UGDF seems reasonable: 

fix, x', kl k',\ t) = y//wz-x.(x',fc;2)/f--(x,fc2,/x2) ■ F{t) , (2.20) 

where fsmau-x{x'-,k!^) is one of the typical small-x UGDFs (see e.g. j26|). So exemplary 
combinations are: KL ® Gauss, BFKL ® Gauss, GBW ® Gauss (for notation see |26| ) . 
The natural choice of the scale is /i^ = M'j^. This relatively low scale is possible with the 
GRV-type of PDF parameterization 27|]. We shall call (I2.20p a "mixed prescription" for 
brevity. 



C. Two-gluon impact factor approach for subasymptotic energies 

The amplitude in the previous section, written in terms of off-diagonal UGDFs, was 
constructed for rather large energies. The smaller the energy the shorter the QCD ladder. 
It is not obvious how to extrapolate the diffractive amplitude down to lower (close-to- 
threshold) energies. Here we present slightly different method which seems more adequate 
at lower energies. 

At not too large energies the amplitude of elastic scattering can be written as amplitude 
for two-gluon exchange [28, 29] 

M,,^,Mt) = ^s^^ I d-h c^.(4)c^.(4) ^^i^"''"|^fj^"'^f^ • (2.21) 

In analogy to dipole-dipole or pion-pion scattering (see e.g. [29]) the impact factor can be 
parameterized as: 

F(kit, k2i) = .0 I n n — \2~^^rVT\ T — \2 ■ (2.22) 

+ (kii + \L2tY + (kij - kat)^ 

At high energy the net four-momentum transfer: t = — (ki^ + kat)^. A in Eg. (12.221) is a 
free parameter which can be adjusted to elastic scattering. For our rough estimate we take 
A = rrip. 

Generalizing, the amplitude for exclusive /o(1500) production can be written as the am- 
plitude for three-gluon exchange shown in FigJD 

^pp-.pp/o(i50o)(s, y, h, t2, 0) = is^^^j^ J (Pkot {as{klt)as{klt))^^^ (as(fc^Jas(A;|J)^^^ 



3F(ko,,ki,)3F(kof,k2,) 



^ fu2 I ..2\fu2 I .,2\fu2 I ,,2^ ^ 99^ M1500) [i^-lt, i<-2t) ■{^■^'i) 
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FIG. 4: The sketch of the two-gluon impact factor approach. Some kinematical variables are shown 
exphcitly. 

At high energy and ?/ ~ the four-momentum transfers can be calculated as: 
tl = -(kot + kli)^ t2 = -(kot - k2tY. 

At low energy and/or y ^ the kinematics is slightly more complicated. Let us define 
effective four-vector transfers: 

92 = {P2-P2) = {q20,q2x,q2y,q2z) ■ (2.24) 

Then ti = qf = qf^ + qf^ and ^2 = ^2= qh + q2f Close to threshold the longitudinal 
components qfi = q^Q — g^^ <^ and = g^Q — <^ 0. Then the amplitude fl2.23p must 
be corrected. Then also four-vectors of exchanged gluons {ko, ki and ^2) cannot be purely 
transverse and longitudinal components must be included as well. To estimate the effect 
we use formula (12.231) ^ but modify the transferred four momenta of gluons entering the 
9*9* ~^ /o(1500) production vertex: 

ki = (0,kit,0) ^ (gio,kit,gi^) , 

k2 = (0,k2t,0) ^ (g2o,k2i,g2.) (2.25) 

and leave k^ purely transverse. This procedure is a bit arbitrary but comparing results 
obtained with formula fl2.23p with that from the formula with modified four-momenta would 
allow to estimate related uncertainties. 

We write the vertex function gg — > /o(1500) in the following tensorial form ^: 

V{ki, k2,PM) = C/o(i5oo)^99 9tivkik2 . (2.26) 



^ It would be more appropriate to calculate in this case a four-dimensional integral instead of the two- 
dimensional one. 

^ In general, another tensorial forms are also possible. This may depend on the structure of the considered 
meson. In principle, the details depend on the form of the vertex. To avoid uncertainties in the kt- 
factorization approach we work in the on-shell approximation. In the on-shell approximation (or infinitely 
heavy meson approximation) the vertex is expressed through decay width and all vertices should be 
equivalent. Even if the off-shell effects are included we do not expect very different energy dependence of 
the cross section for different tensorial forms as due to kinematics only small virtualities of gluons enter 
into game. 
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The normalization factor is obtained from the decay of /o(1500) into two soft gluons: 



|C/o(i50o)^99p - ]vf| (ATS _ i^^ham^gg ■ {221) 
Of course the partial decay width is limited from above: 

r/o(1500)-»gg < ^ toi ■ (2.28) 

The amplitudes discussed here involve transverse momenta in the infra-red region. Then 
a prescription how to extend the perturbative asik^) dependence to a nonperturbative region 
of small gluon virtualities is unavoidable. In the following asik^) is obtained from an analytic 
freezing proposed by Shirkov and Solovtsev 3^ . 



D. Double-diffractive mechanism with intermediate pionic triangle 

The /o(1500) vrvr is the second most probable decay channel 2^. As a consequence the 
mechanism shown in Figj5]may play important role in the exclusive production of /o(1500) 



31 1 . It is relatively easy to estimate the contribution of this mechanism at high energies 
3l|. In this paper we shall make an estimate of the corresponding cross section not far from 



the threshold, where the situation is slightly more complicated. 

Foffik, ki) 

^ /o(1500) 
k ' J D > 

FIG. 5: A sketch of the double-diffractive mechanism with pionic loop for exclusive production of 
the glueball candidate /o(1500). Some kinematical variables are shown explicitly. 

The amplitude of the process pp — »• pp/o(1500) sketched in FigJS] can be written in a 
simplified form [3l| as: 



-MaiA2^a;a!, iy,Plt,P2t,(i)) ~ Tipjpf^{qi,q2,Pfo)5x^X{Fcut{si,eff)5x2X'2Fcutis2. 




(2.29) 

The delta functions are related to helicity conservation in hadronic processes. While the 
pomeron (sub) amplitudes are dominantly imaginary, the reggeon (sub) amplitudes have both 
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real and imaginary parts. The factor rjf ^ i — 1. In the formula above ajp(ti,2) = «jp(0) + 
a'jp ■ ti^2 and ajR(ti,2) = «jr(0) + ■ ti,2 are a so-called reggeon trajectories. For brevity we 
use notation ti^2 which means ti or t2- We take from the phenomenology: ajp(O) = 1.0808, 
a'jp = 0.25 GeV-2, a^lO) = 0.5475 and = 0.93 GeV-^ [3^ The strength parameters for 
the ttN scattering fitted to the corresponding total cross sections [sl]: = 13.63 mb and 
= (27. 56 + 36. 02)/2 mb. * At not too high energies the slope parameter B^^n ~ 6 GeV~^. 
The subchannel Mandelstam variable Si^e// and S2,e// are related to center-of-mass energies 
of relevant pion-nucleon subsystem. In principle, they are functions of pion-four momenta 
in the triangle: Si^e// = ^i,eff{k, ki,p[) and S2,e// = ■S2,e//(fc, ^2,^2) general should be 

put inside of the triangle function Tppj^^^qi, q2,Pfo) which depend on the 3-body kinematics, 
i.e. on four-momenta of the two exchanged pomerons (in general pomeron-reggeon, reggeon- 
pomeron or reggeon- reggeon) . In order to simplify the calculation we take instead average 
values calculated as: 

Sl,eff = {p[+Pfo/2f , 

^2,e// = (P2+P/o/2)'. (2.30) 

The factors Fcut{s) are added to cut off the region of small s, where the naive Regge 
parametrization does not apply. We parametrize them in terms of the smooth function: 

exp ( 

FUS) = \ , (2.31) 

where W = y/I. The parameter Wcut gives the position of the cut and parameter acut 
describes how sharp is the cut off. The latter parameter can have significant influence on 
the numerics. We shall take Wcut = 2 GeV and acut = 0.1 GeV. For large energies Fcut{s) ~ 
1 and close to kinematical threshold W = + Mjy : Fcut{s) ~ 0. This means that we limit 
to double-diffractive contribution only. 

The effective Regge parametrizations of nN interactions [3^ are for both colliding par- 
ticles being on-mass-shell. In our case the triangle pions are off-mass-shell. We correct the 
Regge strength parameters by multiplying by two vertex form factors Foff{k, ki) (see Figj5]). 
We take them in the following factorized form: 

F,;^(fc, h) = exp - ml\/Alfj) ■ exp {-\k^ - mH/Al^j) . (2.32) 

A.off is in principle a free parameter. In the calculation shown in the result section we shall 
take Aoff = 1 GeV. The dependence on triangle four- momenta forces us to merge the form 
factors inside the triangle integration which leads to a modified pion-triangle function: 

/d^k - 
j^T{k; gi, q2,Pfo)Foff{k, ki)Foff{k, ki) , (2.33) 



We take average value for 7r+p and tt p scattering. 
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where standard triangle integrand T{k; qi, q2,Pfo) reads: 



f{k;qi,q2,Pfo) 




1 




^ — ml + ie 

TV ' 

1 



X 



— ml + ie 
1 



X 



F{kl,k2,Pfo) gnnfo 



J^2 _ jj^2 _[_ 



(2.34) 



In addition to three pion propagators we have written three vertex form factors which are 
functions of four momenta of corresponding legs. In principle, these functions are relatively 
well known for space-like pions. We parametrize the triangle-vertex form factors in the 
following factorized exponential form: 



In this factorized form each exponent is associated with individual leg in the vertex. Such 
form factors (exponents) are normalized to unity when pions in the loop are on-mass shell. 
Please note that we symmetrically (modulus in fl2.35p ) damp configurations above and below 
pion-mass shell. A^r is related to the size of the pions in the triangle. It is natural to 
expect: A^^ < Aoff. In the calculation presented here we shall take A,r = 0.5 GeV. Since 
the configurations close to the mass shells give the biggest contributions the sensitivity to 
the actual value of the form factor F (see Eqs. (12.341) and (I2.35P ) is not substantial. The 
gfoTTTT = q-TTTT/n coupllug coustaut can be calculated from the corresponding partial decay 



Calculating the triangle function for running kinematics of the hih2 /i'^/i2/o(1500) 
process (each point of the phase space) is in practice impossible. We calculate numerically 
the triangle function for: 



Transverse components are on average small and are neglected in the present approximation. 
Close to threshold | {qio)y=o \ + \ {.qiz)y=Q \ and | (g2o)y=o I 7^ I {(l2z)y=o \- 

E. Pion-pion MEC amplitude 

It is straightforward to evaluate the pion-pion meson exchange current contribution shown 
in FiglJl If we assume the type coupling of the pion to the nucleon then the Born 



F(gi, ki, k) = exp [~\kl 
F{q2,k2,k) = exp {-\kl 
F{ki,k2,Pfo) = exp {-\kl 



ml\/Al).e^p{-\e-ml\/Al) , 
ml\/Al).e^p{-\e-ml\/Al) , 
ml\/Al).exp{-\kl-ml\/Al) . 



(2.35) 



width 31 . 




(2.36) 
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amplitude squared and averaged over initial and summed over final spin polarizations reads: 



X 



X 



(TP ^^\(F'^^\ ( Pi I 2pi ■ p'l 

r ' + ^ + ^ V(^i + + (^1 + m)(E; + m); J - 




In the formula above m is the mass of the nucleon, Ei, E2 and E[, E2 are energies of initial 
and outgoing nucleons, pi, p2 and p'^, P2 are the corresponding three-momenta and m^r is the 
pion mass. The factor q-^nn is the familiar pion- nucleon coupling constant and is relatively 

■ 2 2 

well known 135] f^^^ = 13.5 - 14.6). In our calculations we take = 13.5. The isospin 



in 



factor Tfc equals 1 for the 7r°7r° fusion and equals 2 for the tt+tt" fusion. Limiting to nucleons 
in the final state, in the case of proton-proton collisions only pp/o(1500) final state channel 
is possible and therefore the vr^vr^ fusion is allowed while in the case of proton-antiproton 
collisions both pp/o(1500) and nn/o(1500) final state channels are possible, i.e. both 7r°7r° 
and 7r+7r~ MEC are allowed. In the case of central heavy meson production rather large 
transferred four-momenta squared ti and ^2 are involved and one has to include extended 
nature of the particles involved in corresponding vertices. This is incorporated via -F7rArAr(ti) 
or -FjrAfAf (^2) vertex form factors. The influence of the t-dependence of the form factors will be 
discussed in the result section. In the meson exchange approach [3^ they are parameterized 
in the monopole form as 



tNN 



it) 



A' 



ml 



A^-t 



(2.38) 



Typical values of the form factor parameter are A = 1.2-1 A GeV i36|, however the Gottfried 
Sum Rule violation prefers smaller A ^ 0.8 GeV 371 . 



The normalization constant |Cp in (12.371) can be calculated from the partial decay width 



as 



\C 



/o(1500)- 



87r 2M|/^„(i5oo)- 



Ml - Ami 



(2.39) 



where rjQ(i5oo)^7r07rO = 0.109 ■ i?i?(/o(1500) — > vrvr) ■ 0.5 GeV. The branching ratio is 
i?-R(/o(1500) 7[tt) = 0.349 [2^. The off-shellness of pions is also included for the 
TTTT /o(1500) transition through the extra Kr7r^/o(i500)(^i; ^2) form factor which we take in 
the factorized form: 



*/o(1500) 



A2 

rnrfo 



ml 



A2 



TTTT Jo 



(2.40) 



It is normalized to unity when both pions are on mass shell 



m,„ 



1 . 



(2.41) 



In the present calculation we shall take A 



1.0 GeV. 
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III. RESULTS 



A. Gluonic QCD mechanism 



Let us start with the QCD mechanism relevant at higher energies. We wish to present 
differential distributions in xp, ti or t2 and relative azimuthal angle 0. In the following 
we shall assume: F 



/o(1500)-^g9 



-ptot 

^ /o(1500) 



This assumption means that our differential 
distributions mean upper limit of the cross section. If the fractional branching ratio is known, 
our results should be multiplied by its value. There are almost no absolutely normalized 



experimental data on exclusive /o(1500) production in the literature, except of Ref . |38 



The absolutely normalized data of the ABCDHW Collaboration [39*] put emphasis rather 
on /2(1270) production. In principle, some (model-dependent) information on glueball wave 
function could be obtained from radiative decays J/ip ^ 7/o(1500) and T — * 7/o(1500) 
(4o| . The present data are not good enough to provide a detailed information on coupling 
of gluons to /o(1500). 

In Figjni we show as example distribution in Feynman xp for Kharzeev-Levin UGDF 
(solid) and the mixed distribution KL (g) Gaussian (dashed) for several values of collision 
energy in the interval ly = 10 - 50 GeV. In general, the higher collision energy the larger 
cross section. With the rise of the initial energy the cross section becomes peaked more and 
more at xp ~ 0. The mixed UGDF produces slightly broader distribution in xp- 




10-^ 



KL UGDF 
KLx Gauss UGDF 

-0.4 -0.2 0.2 0.4 



FIG. 6: The distribution of /o(1500) in Feynman for W = 10, 20, 30, 40, 50 GeV (from bottom 
to top). In this calculation the Kharzeev-Levin UGDF (soUd line) and the mixed distribution KL® 
Gauss (dashed line) were used. 

In Fig 13 we present corresponding distributions in t = ti = t2. The slope depends on 
UGDF used, but for a given UGDF is almost energy independent. 

Finally we present corresponding distributions in relative azimuthal angle between outgo- 
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: fo(1500) : 

1^ 

KL UGDF = 

in-i r 

1^ 0.2 0.4 0.6 0.8 1 

-t (GeV^) 

FIG. 7: Distribution in t = ti = for Kharzeev-Levin UGDF for W = 10, 20, 30, 40, 50 GeV 
(from bottom to top). The notation here is the same as in FigjGl 

ing protons or proton and antiproton ^. These distributions have maximum when outgoing 
nucleons are back-to-back. Again the shape seems to be only weekly energy dependent. 



-e- 

T3 



10^ 
10 



1 

■■Ob 20 40 60 80 100 120 140 160 180 

<[> (deg) 

FIG. 8: Distribution in relative azimuthal angle for Kharzeev-Levin UGDF for W = 10, 20, 30, 
40, 50 GeV (from bottom to top). The notation here is the same as in FiglBl 




fo(1500) 
KLUGDF 



Tlic QCD gluonic mechanism is of course charge independent. 
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B. DifFractive versus pion-pion mechanism 



What about the pion-pion fusion mechanism? Can it dominate over the gluonic mecha- 
nism discussed in the previous subsection? In Figj9]we show the integrated cross section for 
the exclusive /o(1500) elastic production 

pp ^ pfoil500)p (3.1) 

and for double charge exchange reaction 

pp ri/o(1500)n . (3.2) 

The thick solid line represents the pion-pion component calculated with monopole vertex 
form factors 02.381) with A = 0.8 GeV (lower) and A = 1.2 GeV (upper). The difference 
between the lower and upper curves represents uncertainties on the pion-pion component. 
The pion-pion contribution grows quickly from the threshold, takes maximum aX W ^ 6-7 
GeV and then slowly drops with increasing energy. The gluonic contribution calculated 
with unintegrated gluon distributions drops with decreasing energy towards the kinematical 
threshold and seems to be about order of magnitude smaller than the pion-pion component at 
W = 10 GeV. We show the result with Kharzeev-Levin UGDF (dashed line) which includes 
gluon saturation effects relevant for small-x, Khoze-Martin-Ryskin UGDF (dotted line) used 
for the exclusive production of the Higgs boson and the result with the "mixed prescription" 
(KL Gaussian) for different values of the ctq parameter: 0.5 GeV (upper thin solid line), 
1.0 GeV (lower thin solid line). In the latter case results rather strongly depend on the 
value of the smearing parameter. The thick dash-dotted line corresponds to the second 
diffractive mechanism with pionic triangle. It is above the WA102 experimental data point. 
This is probably because of absorption effects not included in the present calculation. This 
contribution stays below the pion-pion fusion contribution at the GSI HESR energies. For 
comparison we show also experimental data point of the WA102 Collaboration from Ref. [3^ 
which lies between the results obtained with "KL" and "mixed" off-diagonal UGDFs. The 
thick long-dashed line corresponds to the second diffractive mechanism with pionic triangle. 
It is above the WA102 experimental data point. This is probably because of absorption 
effects not included in the present calculation. This contribution stays below the pion-pion 
fusion contribution at the GSI HESR energies. 

We calculate the gluonic contribution down to W = 10 GeV. Extrapolating the gluonic 
component to even lower energies in terms of UGDFs seems rather unsure. At lower energies 
the two-gluon impact factor approach seems more relevant. The impact factor approach 
result is even order of magnitude smaller than that calculated in the KMR approach (see 
lowest dash-dotted (red) line in Fig. [9]), so it seems that the diffractive contribution is 
completely negligible at the FAIR energies. 

Our calculation suggests that quite different energy dependence of the cross section may 
be expected in elastic and charge-exchange channels. Experimental studies at FAIR and 
J-PARC could shed more light on the glueball production mechanism. 

The smaller energies the larger values of xi and X2 are involved. Many of unintegrated 
gluon distributions in the literature are formulated in the region of very small x. Extrapola- 
tion of the method down to small energies automatically means going to the region of large 
X. Below we wish to demonstrate this fact. In Fig{TO]we show the ratio of the cross sections 

a{W; Xi < Xo,X2 < xo) 
^= ^) ' ^^-^^ 
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FIG. 9: (Color online.) The integrated cross section as a function of the center of mass energy for 
pp — > pp/o(1500) (left panel) and pp nn/o(1500) (right panel) reactions. The thick solid lines 
are for pion-pion MEC contribution (A = 0.8, 1.2 GeV), the dashed line is for QCD diffractive 
contribution obtained with the Kharzeev-Levin UGDF, the dotted line for the KMR approach and 
the thin solid lines (blue) are for "mixed" UGDF (KL Gaussian) with o-q = 0.5, 1 GeV. The dash- 
dotted line represents the two-gluon impact factor result. The thick long-dashed line corresponds 
to the second diffractive mechanism with intermediate pionic triangle. The experimental data point 
at W = 29.1 GeV is from Ref.|38]. 



as a function of center-of-mass energy. Above xq was introduced to define the region of 
small/large x. The solid line corresponds to a;o = 0.1 and the dashed line to xq = 0.2. Down 
to largest HESR energies one stays in the region of xi, X2 < 0.2. 



C. Predictions for PANDA at HESR 

Let us concentrate now on pp collisions at energies relevant for future experiments at 
HESR at the FAIR facility in GSI. Here the pion-pion MEC (see Figj2]) seems to be the 
dominant mechanism, especially for the charge exchange reaction pp nr;,/o(1500). As 
discussed in the previous section the gluonic component can be there safely neglected. 

In Fig{TT]we show average values of ti (or for the two-pion MEC as a function of the 
center of mass energy. Close to threshold W = 2mjv+wi/o(i5oo) the transferred four-momenta 
squared are the biggest, of the order of about 1.5 GeV^. The bigger energy the smaller the 
transferred four-momenta squared. Therefore experiments close to threshold open a unique 
possibility to study physics of large transferred four-momenta squared at relatively small 
energies. This is a quite new region, which was not studied so far in the literature. 

Below we shall present cross sections for the pp — >■ nn/o(1500) reaction. The cross section 
for the pp pp fo{1500) reaction can be obtained by rescaling by the factor of 1/4. 

The maximal energy planned for HESR is ^/s = 5.5 GeV. At this energy the phase space 
is still very limited. In Fig{T2]we show rapidity distribution of /o(1500) calculated including 
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FIG. 10: The ratio of the cross sections (see Eq. (I3.3p ) as a function of center-of-mass energy. The 
sohd hne corresponds to xq = 0.1 and the dashed line to xq = 0.2. 




10 12 14 16 18 20 
W (GeV) 



FIG. 11: (Color online.) Average value of < ti >=< t2 > as a function of the center-of-mass 
collision energy for the two-pion exchange mechanism. In this calculation A = 0.8 GeV. 



pion-pion fusion only. For comparison the rapidity of incoming antiproton and proton is 
1.74 and -1.74, respectively. This means that in the center-of-mass system the glueball is 
produced at midrapidities, on average between rapidities of outgoing nucleons. 

In FigUHl we show transverse momentum distribution of neutrons or antineutrons pro- 
duced in the reaction pp — > ?T,n/o(1500). The distribution depends on the nNN form factors 
FnNNih) and F^^Nih) in formula fl2.37p . 

In FigHH we show azimuthal angle correlation between outgoing hadrons (in this case 
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y 



FIG. 12: Rapidity distribution of /o(1500) {tt'^'k fusion only) produced in the reaction pp - 
nn/o(1500) for W = 3.5, 4.0, 4.5, 5.0, 5.5 GeV (from bottom to top). In this calculation A = 0. 
GeV. 
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FIG. 13: Transverse momentum distribution of neutrons or antineutrons produced in the reaction 
pp — > nn/o(1500) (tt+tt" fusion only) for W = 3.5, 4.0, 4.5, 5.0, 5.5 GeV (from bottom to top). In 
this calculation A = 0.8 GeV. 

neutron and antineutron) . The preference for back-to-back configurations is caused merely 
by the hmitations of the phase space close to the threshold (the matrix element for pion-pion 
fusion is 0-independent) . This correlation vanishes in the limit of infinite energy. At high 
energy, where the phase space limitations are small, the distributions are isotropic, there 
is no dependence on azimuthal angle. In practice far from the threshold the distribution 
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becomes almost constant in azimuth. This has to be contrasted with similar distributions 
for pomeron-pomeron fusion shown in FiglH] which are clearly peaked for the back-to-back 
configurations. Therefore a deviation from the constant distribution in relative azimuthal 
angle for the highest HESR energy of W = 5.5 GeV for pp — > p/o(1500)p can be a signal of 
the gluon induced processes and/or the presence of subleading reggeon exchanges, e.g. pp. 
It is not well understood what happens with the gluon induced diffractive processes when 
going down to intermediate (W = 5-10 GeV) energies. Our calculations shows, however, that 
the diffractive component is negligible compared to the pion-pion fusion at < 10 GeV. 
Possible future experiments performed at J-PARC could bring some new insights into this 
issue by studying distortions (probably very small) from the pion-pion fusion mechanism. 

-e- 

-D 10^ 



10 



1 



10 b 20 40 60 80 100 120 140 160 180 

4> (deg) 

FIG. 14: Azimuthal angle correlations between neutron and antineutron produced in the reaction 
pp nn/o(1500) (vr+vr^ fusion only) for W = 3.5, 4.0, 4.5, 5.0, 5.5 GeV (from bottom to top). In 
this calculation A = 0.8 GeV. 

Up to now we have neglected interference between pion-pion and pomeron-pomeron con- 
tributions (for the same final channel). This effect may be potentially important when both 
components are of the same order of magnitude. At J-PARC energies there could be, in 
principle, some small interference effect. While the pomeron-pomeron contribution is 
dominantly nucleon helicity preserving the situation for pion-pion fusion is more compli- 
cated. In the latter case we define 4 classes of contributions with respect to the nucleon 
helicities: cc (both helicities conserved), cf (first conserved, second flipped), fc (first flipped, 
second conserved) and // (both helicities flipped). The corresponding ratios of individual 
contributions to the sum of all contributions are shown in FigJTSl In practice, only the cc 
vrvr contribution may potentially interfere with the gluonic one. From the figure one can 
conclude that this can happen only when both transverse momenta of the final nucleons are 
small. We shall leave numerical studies of the interference effect for future investigations, 
when experimental details of such measurements will be better known; but already now one 



At the PANDA energies the problem is rather academic as the diffractive component can be neglected. 
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can expect them to be rather smalL 





FIG. 15: Helicity decomposition of the cross section on the {pit,P2t) plane for W = 10 GeV. i?cc 
(upper left), Rff (upper right), i?c/ (lower left), Rfc (lower right). The standard nucleon dipole 
form factor was used in this calculation. 

Now we wish to show the size of the triangle-double-diffractive (TDD) component at the 
GSI HESR energy range. In Fig{Tn] we compare it with the pion-pion fusion component. 
The TDD component vanishes quickly with decreasing energy and stays below the pion-pion 
fusion component for the HESR energy range. The quick decrease of the cross section is 
caused mainly by the Fcut{si^eff) and Fcut{.S2,ef }) factors in Eq. (12.291) and (12.311) and reflects 
smallness of t^N subchannel energies. 



IV. BACKGROUND FOR THE /o(1500) PRODUCTION FOR pp pp/o(1500) RE- 
ACTION 

A. Pion-pion rescattering background 

In the previous section we have shown that in the PANDA energy range the pion-pion 
fusion is the dominant reaction mechanism for the production of the glueball candidate 
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FIG. 16: Comparison of the pion-pion fusion component (solid lines) and the double-diffractive 
one with pionic triangle (dashed line). The details concerning the double-diffractive component 
are explained in section II-D. The vertical gray band shows the range of the center-of-mass energy 
available by the PANDA experiment. 

/o(1500). Up to now we have calculated the cross section for production of /o(1500) meson 
- a process with three particles (p, p and /o(1500)) in the final state. 




FIG. 17: The 2^4 amplitude for the reaction through the glueball candidate /o(1500). The stars 
attached to vr'^ mesons denote the fact that they are off-mass-shell. 

In practice one must select a given decay channel of /o(1500). There are a few options: 
(a) a two-pion decay (tt+tt^ or tt^tt^), (b) a four-pion decay (c) a two-kaon decay. The first 
one is attractive due to its simplicity but may have a large background. The second requires 
more complicated analysis but may have smaller background. The branching fraction for 
the last option is smaller by a factor of about 5 than for the two-pion channel. 

Let us consider now an estimate of the background to the pp — > ppn^n^ reaction. In 
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FIG. 18: The 2 — > 4 amplitudes with the intermediate p-meson (reggeon) exchange as an example 
of the background. 

FigHZI we present our reaction of interest - the reaction which proceeds through the scalar 
resonance /o(1500). This reaction is viewed now as a process with four particles {p, p, tt"*", tt") 
in the final state. Unavoidably there exists a nonreduceable background to this process. In 
FigHH] we show an example of the background. We shall call these two complex diagrams 
as p- meson (reggeon) exchanges for brevity. The region of Wt^t^ ~ 1.5 GeV is slightly above 
the region of application of the standard meson-exchange formalism and slightly below the 
region of application of high-energy Regge approach. In principle, one should consider both 
approaches. 

In the p meson-exchange formalism the reduced amplitude for the tt^ti^ — > tt+tt^ process 
can be written as: 

,.p-exch. _ p /A (gl + P2.)PfJ.i'{l2 + PI) p /A 

ttOttO— >7r+7r- iJ-mvp-'^ ivKpy'' ) - „ y-KitpJ^ ■K-np\'^ } 

t -mj + iV p-nip 

QTV-irpP-KTrniu] - - ■ — ; 9ttitpFtyttp(u) . (4.1) 

u — mj + iTpUip 



Above 



F^^p{k') = exp ( — ^ ) = exp ( ^{k' - m^) ) . (4.3) 



Pi.u{^) = -Qf^v + kpkjml . (4.2) 

The quantities F.„T,p{k'^) in (14.11) describe couplings of extended objects: pions and the ex- 
changed p-meson. We parameterize them in the exponential form: 

Consistent with the definition of the coupling constant the form factors are normalized to 

unity when p meson is on-mass-shell. We take = 2.6, which reproduces the p meson 
decay width [2% and A = 1 GeV (E^^p = 4 GeV"^). 

In the case of p-reggeon exchange the amplitude can be written as 

MlZo'ZZ- = S3,Vp{i)C.npF{i) F{i) + suVp{u)C^npF{u) ( ^ j F(tx) . 

(4.4) 
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We parameterize the vertex form factors in the standard exponential form used usually in 
the Regge phenomenology 

F(e) = exp(^kA . (4.5) 



In practical calculation we take -Bttttp = 6 GeV^^. In the formula (14.41) rip{k^) is a signa- 
ture factor which we take here ?7p ~ i + 1 and ap(t/u) = «p(0) + a'^ ■ t/u is a so-called 
reggeon trajectory. We take from the phenomenology: ap{0) = 0.5475 |32| and a'p = 0.9 
GeV~^. The strength parameter C^rTrp can be obtained assuming Regge factorization (see 
e.g. [2^) and using the known strength parameters for the A^A^ and vrA^ scattering fitted to 
the corresponding total cross sections [32]. The simple Regge parameterizations apply for 
energies W > 2 GeV (see e.g.js^]). In our case of pion-pion scattering energies of W ^ 1.5 
GeV are of interest. Here a small modification of the Regge formula (14. 4p may be in order. 
Consistent with meson-exchange formalism (spin-1 exchange) one may expect saturation of 
the tt^tt" —>■ 7r"'"7r~ cross section at lower energies. The following freezing of the energy factor 
in (14.41) seems a reasonable correction: 



So J \ So 

t2 



(4.6) 



where Sfreez — Wf^^^^. One may expect Wfreez = 1-5 - 2.0 GeV. The compatibility of 
the Regge formalism with low-energy approaches for pion-pion scattering was discussed in 



Ref. 33 



The 2 — i> 2 amplitudes (14.11) and (14. 4p may be inserted into the 2 — >^ 4 amplitude of Fig{ 
When doing so we include in addition the correction for off-shellness of incoming pions as 
was done for the /o(1500) meson using exponential form factors (as in 14. 9p . Now we can 
perform a genuine 2 — > 4 calculation including four-body phase space. 

Here we discuss the results with p-reggeon exchange only. We have checked that the p- 
meson exchange formalism discussed in this section provides approximately the same results 
at Wtttt < 1.2 GeV as the modified reggeon exchange with Wfreez = 1-5 - 2.0 GeV (see 
formula (I4.6p ). Therefore the modified reggeon-exchange calculation provides a realistic 
predictions in the broad range of pion-pion energies, both above and below the /o(1500) 
resonance. 

In Fig {19] we show two-pion invariant mass distribution. The solid line corresponds to our 
resonance contribution. The dashed lines correspond to the p-reggeon exchange contribution. 
Here the resonance contribution is much lower than the p-exchange background. In this 
calculation the integration over whole phase space was done. 

In Figj20]we show corresponding distributions in rapidity of vr"'" or vr" (identical). In order 
to better see the overlap of the signal and background for the p-reggeon exchange we impose 
in addition 1.4 GeV < M^-^ < 1.6 GeV (the region of the /o(1500) resonance). Limiting to 
very small center-of-mass rapidities one can further improve the signal-to-background ratio. 
In Figj2T]we show two-pion invariant mass distribution with extra cuts: -0.5 < U-k+iU-k- < 
0.5. While the /o(1500) contribution is only slightly modified, the p-reggeon background 
contribution is reduced by more than order of magnitude. One can clearly see the signal 
over background in this case. Especially the high-energy side of the /o(1500) meson is now 
free of the p-exchange background. A better separation can be done by using pion-pion 
partial wave analysis. 

Close to the two-pions production threshold the Roper resonance excitation and its sub- 
sequent decay (A^*(1440) — > Ntttt) is known to give the dominant contribution to the 
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FIG. 19: Pion-pion invariant mass distribution for /o(1500) (solid line) and the p-reggeon exchange 
background (dashed lines) for naive (upper line) and corrected (two lower lines with Wfreez = 1-5, 2 
GeV) extrapolations to low energies. Here the full phase space has been included. The calculation 
was performed for the highest PANDA center-of-mass energy W = 5.5 GeV. 




FIG. 20: Rapidity distribution of pions from the decay of the glueball candidate /o(1500) (solid line) 
and from the p-reggeon exchange background (dashed lines) for naive (upper line) and corrected 
(lower lines) extrapolations to low energies. For the background we impose in addition: 1.4 GeV 
< Mtttt < 1.6 GeV. The calculation was performed for the highest PANDA center-of-mass energy 
W = 5.5 GeV. 
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FIG. 21: Pion-pion invariant mass distribution for /o(1500) (solid line) and the /?-reggeon exchange 
background (dashed lines) for naive (upper line) ajnd. corrected (lower lines with 'l^freez — 

1.5, 2.0 

GeV) extrapolations to low energies. Here an additional condition on center-of-mass rapidities 
-0.5 < ?/7r+,y7r- < 0.5 has been imposed. The calculation was performed for the highest PANDA 
center-of-mass energy W = 5.5 GeV. 



pp ppiT^iT reaction S^l- The same may be expected also for the pp ppn^Ti reaction. 



The Roper resonance produces the two-pions in dominantly the 1=0 and /=0 state (the 
tail of the a meson), i.e. the strength is concentrated at Mt^^ much lower than /o(1500). The 
kinematical constraint gives M^^^ < M^r* (1440) — Mj^ ^ 0.5 GeV. In addition, this contribu- 
tion could be eliminated by extra cuts on invariant masses M{p7c^TT~) and M{pTT~^7i~). The 
same method can, at least in principle, be used to eliminate the double A, A excitations 
followed by their decays A ^ irp and A ^ vrp In this sense the last two contributions 
(Roper and double isobar excitations) are reduceable. To which extend precision of the real 
apparatus will allow such a reduction is a matter of further investigations. 

Certainly complete analysis requires including more processes and an analysis of cuts 
allowing for improving the signal-to-background ratio. This certainly goes beyond the scope 
of the present paper where we only signal a huge increase of the /o(1500) meson production 
cross section in the PANDA energy range due to pion-pion fusion, the mechanism never 
discussed before in the literature. 

The pp — > PP71717171 reaction may be more favorable as far as the signal-to-background 
ratio is considered. Unfortunately theoretical calculation of background are not feasible in 
this case. It is not clear to us at present if the 6-body channel can be measured by the 
PANDA detector at FAIR. 



Here I is angular momentum between pions and / is the total isospin of the pion pair. 

We have checked that eliminating the region of double-Delta excitation at the highest PANDA energy W 

= 5.5 GeV reduces the signal by less than about 5%. 
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FIG. 22: A sketch of the double-diffractive mechanisms of exclusive production of the vr+vr pairs. 



In the case of the pp — > ppK^K~ reaction the relevant branching fraction is smaller, 
but not negligible (/o(1500) KK = 8.6 % [23|). On the other hand the contribution 
from nucleon resonances is probably considerably smaller. There is, however, unreduceable 
contribution from K* exchange in the K^K^ K^K~ subprocess. The parameters for the 
latter reaction are much less known than those for the tt^tt^ tt+tt^ subprocess. 



B. Double-difFractive two-pion background 

At high energies another two-pion continuum may be of interest. The underlying mech- 
anism was proposed long ago in Ref. 4l|]. The general situation is sketched in Fig 1221 The 
corresponding amplitude for the pp —>■ pprc'^Ti^ process (with four-momenta Pa + Pt 
Pi + P2 + P3 + Pa) can be written as 

+ M^^iti, su) F{tb) ^ ^ , Fih) M23(t2, S23) , (4.7) 

where Mj^ denotes "interaction" between nucleon i=l (forward nucleon) or i=2 (backward 
nucleon) and one of the two pions k=3 (vr"*"), k=A {it'). In the Regge phenomenology they 
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can be written as: 



(/„ \ QjR(tl)-l / N Qjp(tl)-1 



\ / \ Qip(tl)-1 

M,4(ti, .u) = ^.14 I ( *^ + Cjp ( 



M,3(t„ .,3) = ^^23 I C'i (^) + Cjp (^j M . (4.8) 

Above Sik = W^i^, where Wik is the center-of-mass energy in the {i, k) subsystem. The first 
terms describe the subleading reggeon exchanges while the second terms describe exchange 
of the leading (pomeron) trajectory. We have neglected real parts of the reggeon exchanges 
amplitudes for simplicity. The strength parameters of the vrA^ interaction are taken from 
Ref.Q. 

The extra form factors F{ta) and F{ti,) "correct" for off-shellness of the intermediate pion 
in the middle of the diagrams shown in FigJ2^ In the fallowing they are parametrized as 



t — ml 



m = exp I 1 , (4.9) 



"off 



i.e. normalized to unity on the pion- mass-shell. We take Aoff = 1 GeV. More details of the 



calculation will be presented elsewhere |42| . The 2 — > 4 amplitude (14. 7p is used to calculate 
the corresponding cross section including limitations of the four-body phase-space. 

To excludes resonance regions we "correct" the Regge parametrization (14.81) by multiply- 
ing byfactors Fcut{sik) (as in 12.311) . In Fig 123] we show the two-pion invariant mass distri- 
bution of the double- diffract ive background together with the /o(1500) signal (vrvr fusion). 
We show three curves corresponding to different cuts on both W^n- Wmin = 2.0 (solid), 
1.9 (dashed), 1.8 (dotted) GeV. The figure suggests that the double-diffractive background 
should not disturb observing the /o(1500) signal at the PANDA experiment energies. 

In this case, unlike for the two-pion involved p exchange discussed in the previous sub- 
section, imposing cuts on pion rapidities would not be helpful as the double-diffractive 
contribution is concentrated at midrapidities as shown in FigJ2H At lower HESR energies 
the situation is better, the double-diffractive two-pion background is relatively smaller. 



V. DISCUSSION AND CONCLUSIONS 

For the first time in the literature we have estimated the cross section for exclusive 
/o(1500) meson (glueball candidate) production not far from the threshold. We have included 
both gluon induced diffractive and triangle-double-diffractive mechanisms as well as the 
pion-pion exchange contributions. 

The QCD diffractive component was obtained by extrapolating down the cross section in 
the Khoze-Martin-Ryskin approach with unintegrated gluon distributions from the literature 
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FIG. 23: Pion-pion invariant mass distribution for /o(1500) (solid line) and the double-diffractive 
background (dashed lines) for different (sharp) cut-off for W-j^jy. Here the full phase space has been 
included. The calculation was performed for the highest PANDA center-of-mass energy W = 5.5 
GeV. 




FIG. 24: Rapidity distribution of pions from the decay of the glueball candidate /o(1500) (solid 
line) and from the double-diffractive background (dashed lines) for different (sharp) cut-off for 
WttN- For the background we impose in addition: 1.4 GeV < M^j^ < 1.6 GeV. The calculation was 
performed for the highest PANDA center-of-mass energy W = 5.5 GeV. 

as well as using two-gluon impact factor approach. A rather large uncertainties are associated 
with the QCD diffractive component. At present only upper limit can be obtained for the 
diffractive component as the /o(1500) — » gg decay coupling constant remains unknown. The 
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coupling constant could be extracted only in high-energy exclusive production of /o(1500) 
where other mechanisms are negligible. 

We have found rather large contribution of pionic-triangle-double-diffractive component 
at higher energies {W > 10 GeV). However, at the GSI HESR energies this contribution is 
strongly damped because of the phase space limitations on the ttN subchannel energies. 

A future experiment at RHIC could contribute to shed some light on the competition of 
the both diffr active mechanisms. 

The calculation of the MEG contribution requires introducing extra vertex form factors. 
At largest PANDA energies they are relatively well known and the pion-pion fusion can be 
reliably calculated. The situation becomes more complicated very close to the threshold 
where rather large and \t2\ are involved. The cross section for energies close to the 
threshold is very sensitive to the functional form and parameters of vertex form factor. 
Therefore a measurement of /o(1500) close to its production threshold could hmit the so- 
called ttNN form factors in the region of exchanged four-momenta never tested before. 

We predict the dominance of the pion-pion contribution close to the threshold. Our calcu- 
lation shows that the diffractive components (in fact its upper limit for the QCD mechanism) 
are by more than order of magnitude smaller than the pion-pion fusion component in the 
energy region of future PANDA experiment. 

The diffractive components may dominate over the pion-pion component only for center- 
of-mass energies W > 15 GeV. Taking into account rather large uncertainties the predictions 
of this component should be taken with some grain of salt. Clearly an experimental program 
is required to disentangle the reaction mechanism at energies 1^ > 15 GeV. 

Disentangling the mechanism of the exclusive /o(1500) production not far from the me- 
son production threshold would require study of the pp pp/o(1500), pp nn/o(1500) 
processes with the PANDA detector at FAIR and pp — > pp/o(1500) reaction at J- 
PARG. In the case the pion exchange mechanism is a dominant process one expects: 
(t{pp — > nn/o(1500)) = 4 x a{pp — > pp/o(1500)). At high energies, when the gluonic, or 
double-diffractive with intermediate triangle, components dominate over MEG components 
a{pp pp/o(1500)) > a{pp nn/o(1500)). 

At intermediate energies one cannot exclude a priori sublcading reggeon exchanges like pp 
for instance. However, we do not know how to reliably calculate them from first principles. 
We believe that the distortions from the pion-pion at low energies and/or distortions from the 
QCD gluonic mechanism at high energy may tell us more and allow for a phenomenological 
analysis taking into account the pp component explicitly. We leave this problem for a future 
analysis when experimental data will be available. 

Only a careful studies of different final channels in the broad range of energies could help 
to shed light on coupling of (nonperturbative) gluons to /o(1500) and therefore would give a 
new hint on its nature. The experimental studies of exclusive production of /o(1500) are not 
easy at all as in the tttt decay channel one expects a large continuum. We have performed 
an involved calculation of the four-body ppTr'^ir' background. Our calculation shows that 
imposing extra cuts should allow to extract the signal of the glueball /o(1500) candidate at 
the highest PANDA energy. A partial wave tttt analysis should be helpful in this context. 
The two-pion continuum will be studied in more detail in our future work. A smaller 
continuum may be expected in the KK or four-pion /o(1500) decay channel. This requires, 
however, a good geometrical (full solid angle) coverage and high registration efficiencies. 
PANDA detector seems to fulfill these requirements, but planning real experiment requires 
a dedicated Monte Carlo simulation of the apparatus. 
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It is a central problem of our field if /o(1500) is a qq or glueball type. Unfortunately, 
our analysis does not allow to give a definite answer to this important question. Some 
information on baryon-baryon correlation may be helpful but certainly not decisive. 

If the cross section at high energies (where the contribution of subleading reggeon ex- 
changes may be neglected) is much smaller than predicted based on the KMR method it 
means that gluons only weakly couple to /o(1500). This could provide some indirect informa- 
tion on the /o(1500) structure. A direct comparison of the shape of differential distributions 
at high energies may provide a valuable test of the KMR method originally proposed for 
exclusive Higgs production (the latter experiment is very difficult as very small statistics is 
predicted). A possible disagreement with the prediction for exclusive /o(1500) production 
at high energies could put into question the KMR approach, at present state of art in the 
field. Experiments at RHIC could be useful in this context and could shed light on the 
nonperturbative couphng of gluons to /o(1500). 
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